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FOREWORD 

This document is one of two documents that constitute the 
final report for Contract NAS8-28089, "Study of Viscous Mixing 
Plume Flow Field." This study was performed by the Lockheed- 
Huntsville Research & Engineering Center, Inc., for the National 
Aeronautics & Space Administration, George C. Marshall Space 
Flight Center, Alabama. 

The other document is "Mathematical Representations of 
Turbulent Mixing," LMSC-HREC D306102, January 1973. 

The NASA Contracting Officer' s Representatives (COR) for 
this contract are Dr. T. F. Greenwood and Mr. D. C. Seymour, 
S&E-AERO-AT. 

The authors extend appreciation to J. A. L. Thomson, 

W. F. Herget and in particular, J. E. Reardon for the helpful 
suggestions and comments. 
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Section 1 

INTRODUCTION AND SUMMARY 

A spectroscopic instrumentation system was developed by Rocketdyne 
Division of North American Rockwell which was used to measure tempera- 
ture and concentration distributions in (hopefully) axisymmetric and two- 
dimensional combusting flows (Ref. 1). This measurement technique has 
become known as zone radiometry. 

The success of the method depends both on how accurately the detected 
radiation can be converted by analysis into the desired temperatures and 
concentrations and on how closely the flow meets the dimensional limitations 
of this measurement scheme. Since the technology of radiative transfer was 
being very actively researched during the same time that the zone radiometry 
experiments were being performed, a critique of the Rocketdyne data redcution 
procedures is in order to determine whether or not application of the present 
state-of-the-art radiation analyses can yield more accurate flowfield 
information. 

Theoretically, a temperature and a partial pressure distribution for a 
given species can be determined from a set of measurements made at one 
particular spectral level. If sets of measurements at more than one spectral 
level are made, partial pressures of several species and some average tem- 
perature can be determined. Practically, such multiple determinations have 
not yet been possible. At best, one set of measurements has been used to 
establish temperature and the partial pressure of one species, and another 
set used to determine the partial pressure of a second species. Therefore, 
an analysis at one particular wavelength is all that is required of available 
zone radiometry data. Such an analysis is described in this report. 

The final goal of this report is to present a recommended data reduction 
scheme for the zone radiometry system. The limitations in this scheme will 
be clearly stated and quantitatively evaluated when possible. 
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To appreciate the utility of zone radiometry methods, one should realize 
that the technique was developed and used extensively to measure "axisym- 
metric" rocket motor plumes. All propellant systems cannot be measured 
with satisfactory accuracy by this method. Plumes with carbon particles may 
become too "optically thick" for the transmission part of the measurement to 
be made, and plumes with only water vapor as the optically active species may 
be too "thin" for accurate measurements. Furthermore, all motors lack axial 
symmetry to some degree; no good measure of this feature has yet been devised. 

Despite the fact that much experimental data are available, no definitive 
comparison with calculated flows exists. Experiments with carefully designed 
burners are the primary source of radiation property data, but when a real 
rocket motor is studied both analysis and experiment become more difficult. 

Zone radiometry has also been used on two-dimensional mixing studies. 
These studies have not been as extensive as those on motors. Again definitive 
comparisons with calculations have not been made, nor have error analysis 
for reducing data been previously reported. 

This report serves as a prelude to a more complete data comparison 
study which will be forthcoming. A detailed treatment of the radiation analysis 
and synopsis of the zone radiometry method is reported herein, so that questions 
regarding the accuracy of reported data and experiments can be determined. 
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Section 2 

RADIATION ANALYSIS 

The zone radiometry system is an instrument used to determine temper- 
ature and composition of optically active species in either an axisymmetric or 
planar gas flow field. The radiative energy transfer analysis which must be 
used to relate the radiation measurements to the flow properties is described 
in this section. 

The basic radiative exchange process is sketched in Fig. 1. Radiation 
from the hot zones is seen by the detector at all times. When the chopper 
is open, radiation from the source which is even hotter than the zones plus 
radiation from the hot zones is seen. Particular frequency (or wave length) 
intervals are measured with the detector by selectively filtering away the 
unwanted radiation. The radiative exchange process is essentially one- 
dimensional, as the angle (3 is quite small; hence, one line of sight is viewed. 
The dimensions of the gases on the immediate sides of 3 are assumed to be 
such that radiative equilibrium with the adjacent, lateral-gas zones is main- 
tained, so there is no net radiative exchange in the lateral direction. 

To quantitatively describe radiative exchange, the concept of intensity 
must be used. At a point, P, consider the monochromatic intensity, 


I = Limit 

dcr, dft, dv 
dt — ► 0 



( 2 . 1 ) 


The term represents the radiant energy in [v, v + dv), where v is the 
frequency of the radiation. The term t is time; the geometric factors cr, 
0, are defined by Fig. 2. The fact that the above limit exists is an 
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Chopper 


Fig. 1 - Zone Radiometer System 



dcr = incremental surface 

N = surface normal 

N1 = parallel to surface 
normal 

d$7 = increment of solid 
angle 

0 = angle between solid 

angle direction and 
surface normal 

(3 = solid angle of all dT2 

ove r da 

(da is located at point R) 


Fig. Z - Definition of Geometric Terms Used to Define Intensity 


4 

LOCKHEED - HUNTSVILLE RESEARCH & ENGINEERING CENTER 



LMSC-HREC D306101 


experimentally observed fact (Milne, Ref. 2, p. 84). Radiation is emitted 
from each point on do; therefore, an integration in is required to calculate 
the radiation flux through da. In optics, this is not the case as intensity is 
defined at a point with da missing in the limit expression. 

1^ is independent of S unless it is modified by the transmitting medium, 
whereas E^ is not. There are other intensities which could have been defined; 
they are I and In . These are defined on the basis of a unit of wave number, 

0), or wave length, X, in the limiting expression. If the transporting medium 
has a unit or known index of refraction, conversions between these intensities 
can be easily made. I^is somewhat more basic because it is independent of 
the index of refraction. However, the overriding criteria to use in selecting 
the intensity to use is the availability of property data. These data are avail- 
able in select wave number increments; therefore, I will be used. 

Radiation in the absence of emission is attenuated according to 


X co l S ii = ^i^i) ex P 


i+l 

/ 


K coP dS 


( 2 . 2 ) 


where here and henceforth brackets indicate functionality and where is 
mass absorption coefficient and p is the density of the absorbing medium. 

If and p are independent of S this relationship is called the Beer-Lambert 

®2 

law. In general, optical thickness = j* K^pdS. A spectral absorption coef- 

S 1 

ficient may be defined by the equation shown on the following page. 
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_ \o^i+l^ *00 _ h* (absorbed) 

W *00 ^i+ 1 ^ *00 (incident) 


(2.3) 


In general, radiation may be absorbed, reflected or transmitted; or, fractionally, 


a + p + T =1 

0) ^00 CO 


(2.4) 


Thus 


V s i> 

r = 1 - a = - T.+ t~ / q ' \ , ifftsO 

w 03 *oo ( ®i+l w 


(2.5) 


or 


T co = exp (/" K 0) P dS ) 


( 2 . 6 ) 


The two absorption coefficients are related by: 


a 0) = 1 - ex P (/ ■ K 0)P dS ) 


(2.7) 


Elements along the solid angle will not only absorb radiation but will 
emit at a rate of 


dE 


00 


dt 


= (p dcr dS) dco df2 


( 2 . 8 ) 


where J 4 is the emission coefficient. J . will be isotropic. 
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The geometry is such that the detector is normal to the view angle 
through the plume; hence, cos 9=1. 

2.1 THE EQUATION OF TRANSFER 

Now a radiation heat balance on a control volume consisting of the solid 
angle P between and ^ can be made. Consider three cross sections of 
P, those at S^, S^^, Sj+j/ 2* call them A^, A.^_j, ^i+l/2 ' Since we wish to 
calculate the radiation to the detector, let the radiation at ^ be I^and that 

at S. be I /v + d I , i.e. I is positive in the negative S direction. The heat 

1 GO 00 00 r ° 

balance becomes: 


I + dl dn A., . dt = 
0 ) 0 ) / 1+1 


A,* 1 A i dt -/( I co +dI o,) 

2 Q 

- / K U ' 5 ( S i + l- S i) I co dndtA i + l/2 

£2 

+ / J M P A i + l/ 2 ( S i + l 


(2.9) 


The limits on the Q integration are over the solid angle P ; all variables are 
constant with respect to this integration. The integration converts intensity 


to flux (Milne, Ref. 2, p. 85). Since P is small and S. 


S. + 1 , Eq. (2.9) 


becomes, 


I A 

00 


p p 

i dt /dSJ -(l^+dl^A^dt J 6a = 


- K p dS I dt A. 

60 r 


00 


L i+l/2 




<K2 


+ P A i+ 1/2 38 


dt 


/ 


d£2 


( 2 . 10 ). 
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Let 

A i+l/2 = A i + T 

\ ( A i + l/2 - J t) * * (h* + dI M )( A i + l/2 + J r) dt P 

= ' K W P dS L dt A i+l/2 P 

+ J 6)P A i+ 1/2 dS dt P 

Dividing by dt (3 dS and neglecting products of differentials. 

Assume the flow to be in local thermodynamic equilibrium, such that 
theory can be used to give 

J co ^"00 ^ cob 

where is Planck’s blackbody intensity 

2 C 2 *o) 3 

A cob f , fi Cco . ,1 

l ex p ( “sr ) ' 1 J 

Similar intensities based on other measures of spectral interval may 
defined. For convenience, several of these are tabulated in Table 1. 

dx = p dS 
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(2.13) 
Kirchoff's 

(2.14) 

(2.15) 

also be 
Let 

( 2 . 16 ) 
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Table 1 

RADIATION RELATIONSHIPS WITH RESPECT TO FREQUENCY, 
WAVELENGTH AND WAVE NUMBER 


I. INTENSITIES (SIEGEL AND HOWELL, REF. 3, PP. 20 AND 31 - FOR 
THE INDEX OF REFRACTION EQUAL TO ONE) 


2 £ v~ 


c 2 (i 

[ ex p (ff)] " *) 

7T 

1 


2/b C 2 

N Xb 

> 5 (l 

[ exp (* at)] ' 1 

)' * 


Z 3 
2 h C CO 

. N ub 


([«*(¥*)] -•) 


where V is frequency in (time ), X is wavelength in (length), and CO is wave 
number in (length *). C is the speed of light in a vacuum, ft is Planck's 
constant, and ^ is Boltzmann's constant. 


(Continued) 
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Table 1 - (Continued) 


II. 

ILLUSTRATION OF v, X AND 00 RELATIONSHIPS (SIEGEL AND HOWELL, 
REF. 3, p. 20) 


For Carbon Dioxide ( CO 2 ) 


For Water (H 2 0) 

Xv 

= C (Speed of light in a vacuum) 


\v = C (Speed of light in a vacuum) 

for 

X Q = 4.45/lx for CO z 

for 

X Q = 2.49 JU. for H 2 0 


V - — - = 6.736 x 10 12 sec 1 

n a 


V = -p- = 1.2039 x 10 14 sec 1 


v A 

O 


° A o 


it 

0 

3 


X 0) = 1 

o o 


00 = -r— = 2247 cm 1 

n A ■ 


oo = = 4016 cm 1 


o 


° o 

f° r 

doo = +25 cm ^ 

for 

doo = +25 cm 1 


dX = doo = -1.98 x 10~ 7 doo 

r, o 


1 _ Q 

dX = 5 " doo = -6.20 x 10 doo 

00 

o 


dX = - 0.465 x 10~ 5 n 


dX = - 1.550 x 1 0 - 6 /lx 


16 

dv = - dX = -1.514 x 10 

dX 

V o 17 

dv = - — - dX = -4.834 x 10 w dX 


dv = +0.705 x 10 +11 sec' 1 


dv = + 0.748 x 10 +12 sec' 1 
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Then the energy balance becomes 


d I 


CD 


K dx 

CO 


X co ^cob 


(2.17) 


Equation (2.17) is called the equation of transfer and is of fundamental 
importance in radiative transfer. The derivation given is consistent with 
Milne (Ref. 2). Other discussions of this equation are given by Viskanta 
(Ref. 4) and Kourganoff (Ref. 5). 


The detector in the zone radiometer system does not indicate I , but 

J 00 

rather the product (0 1^). Since 0 is a constant, the signal is proportional 
to 1^ . 0 is a definite number, the view angle of the radiometer. The im- 

portant point is that Eq. (2.17) is valid for any constant value of 0 . More 
will be said of these choices in subsequent pages. 


To solve the equation of transfer, an integrating factor is introduced 
so that the two terms on the LHS of Eq. (2.17) may be combined. 




v 1 ' ( -/ 




* hub exp 


A. 

f- 




(2.18) 


Primes denote dummy variables. 

x 


V cx pI / - 




.o 

f 


x CDb 

o \ o 


f 


L,K<=*P( / - K U <fa ") K U dX ' 


(2.19) 


= U x I exp ( / - K u *=) + / exp ( / 

'o ' O \ o 


/ • K co dx " ) K co dx ' 


(2.20) 
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This equation affirms that the intensity which arrives at 0 comes from x 
and is attenuated between x and 0 or from emission plus self-absorption 
along x to 0, thus the two terms on the RHS of Eq. (2.20). 

Note the first term on the RHS of Eq. (2.20) is often omitted with the 
understanding that x becomes so large that the path becomes '’optically 
thick.” This point is discussed in Goody (Ref. 6) p.456. When this omission 
is used, the term is recovered by using a boundary condition on the RHS, 
when the integration is performed, that equals the omitted term. Such a 
ploy will be used here. 

1^ is a monochromatic radiation intensity. Experimentally, a specific 
wave number cannot be isolated for study, so a wave number interval is used. 
An appropriately averaged intensity, called radiance , is obtained by 



(2.23) 


(2.24) 


Therefore 



= -h f f - "ab (~5r) dx ' d “ 

col O 
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To exchange the order of integration, define 


and 


co2 


= -L f 
aco y 

col 


T dCO 
CO 


(2.25) 


0)2 


'cob 






(2.26) 


col 


Then 


CO 


/ 

/ 




(2.27) 


To be consistent r should have a subscript u, but this is not convenient. 

Golden (Ref. 7) strongly contested the possibility of this inversion of 
order. Simmons (Ref. 8) presented several other ways of accomplishing the 
averaging and an alternate derivation. 

Equation (2.27) may be approximated with finite -differences as: 


N 

£ 

i = 1 




h} [vi - 


T i 


(2.28) 


or 


V I oobl (1 “ T l ) + X cob2 (T 1 " T 2 ) + * * • ^obN^N-l ' T N ) 


(2.29) 


since 


a =0, t = 1. 
o ’ o 
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It is crucial to the understanding of non-grey, radiation analysis to 
appreciate that the quantities in Eqs. (2.27), (2.28) and (2.29) are spectrally 
averaged; however, from this point on the overbars will not be shown because 
only averaged quantities are of interest. 

2.2 THE EVALUATION OF T AVERAGE 

Equation (2.28) represents the intensity which is proportional to that 
detected in the zone radiometery experiments. The summation can be directly 
evaluated, but even though the blackbody intensities are well behaved the trans- 
mittances are not. This is the reason that Krakow et al. (Ref. 9) developed the 
Curtis- Godson approximation to properly average the t's. The development 
of this approximation is given below. 

Two problems must be simultaneously solved. First, an average over 
a certain wave number interval must be established, because radiative trans- 
itions occur in a discontinuous manner even with respect to the narrow accep- 
tance interval of the detecting device. For isothermal, homogenous gases, 
this wave number averaging has been accomplished with "band models." The 
Random Band model with constant line widths was chosen for this purpose. 

It is stated: 

- In r = 2 tt f |X[ (2.30) 


where { } denotes functionality and 



(2.31) 


and are band model parameters and f is the probability distribution 

for line strengths. S is still distance. Specific values of the band model 
parameters will be subsequently quoted. 
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The second problem which must be overcome is to devise a way that 
the band model representation of the transmittance may be used for a n on- 
homogeneous, nonisothermal path. The Curtis -Godson approximation may 
be used to calculate an average for such paths. Let 

-(f) s - £ (f ). s h < 2 - 32 > 

7 h h 

and 

(fMf)- £ (f) h s h (i) h < 2 - 33 > 

h 

where h represents a zonal increment of constant temperature and composition. 
By combining Eqs. (2.31) and (2.33) 

x = [L (s/d) h S J / 5jr Z (s/d) h s h ( ^ /d) h (2 34) 

1 h J ' h 


Within each zone 


X* = (s / d) h S h^ Zn (T / d) h (2,35) 

- In T* = 2tt (y/d) h f (2.36) 


where stars emphasize zonal properties. 
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Therefore, Eq. (2.30) becomes 


-in r 


£ < (-In T* h /' {x*}) 

E< (-*" </ f {<}) ’ 


[E <(-«“</« {<})] 

E <(-«»</' K }) 2 ' 


Equation (2.37) was derived and substantiated by experiments in Krakow et al. , 
(Ref. 9). Two limits for this expression exist. 

If X — ► 0, i.e., is less than 0.1, 


If X 


-£n r = 2 ( _in T h)‘ 

h 

oo, i.e., is greater than 3, 


-£n t = 


EKO 2 


l/2 


(2.38) 


(2.39) 


Now the equations developed in the two previous paragraphs may be used, 
if appropriate band model data are available. The General Dynamics experi- 
ments (Refs. 10, 11 and 12) provide such data. An exponential probability distri- 
bution was used and values of (s/d) and (y/d) were determined. 


f|x| = X 


1 + 


n 



- 1/2 


(2.40) 
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~fn T = 


(s/d) S 


’i+JaZdl s 

4 (y/d) 


TFT 


(2.41) 


Furthermore, let (s/d) = k and (y/d) = a 


-j£nr = 


kS 


23 k h S 1 

h 


1 + 


kS 

4a 


1/2 



(? k h\) 2 1 

1 + - h 

1/2 


I T 

4 2 k h S h a h 



(2.42) 


This is the same relationship that is used by Reardon and Huffaker (Ref. 13) 
to calculate radiation from a line of sight. For a single isothermal, constant- 
composition zone: 


* 

-In t 


k S" 


1 + 


kS 

4a 


772 


(2.43) 



(2.44) 


P Q = reference state of 1 atmosphere 
k is in (cm *) and is tabulated for H 9 0, CO and 

O Ld 

CO 2 in the General Dynamics reports (Refs. 10, 11 and 12). The term 
represents the partial pressure of the radiating species in atmospheres. 
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Unless the pressure is much lower than one atmosphere, Doppler 
broadening is negligible with respect to collision broadening. Assuming 
such a case, "a" can be calculated from the tabulated data in Ref. 12, pp. 22 
23 or from Reardon and Huffaker (Ref. 13) pp. 141-144. 

For continuous radiators, grey gases (throughout the spectral range 
of interest), "a" — oo and 


-£n T = k S (2.45) 

This correctly implies that the optically thin limit Eg, (2.38), can be used to 
calculate integrated values of r . Thus non-grey gases, which are optically 
thin because of geometry and density distributions in addition to grey gases 
obey Eq. (2.38). 

In general, "a" is the. fine structure parameter which is the ratio of 
the line width, y, to line spacing, d. Width of a radiating line is broadened 
by collisions between the atoms and molecules of the gas. The general form 
of the line width with collisional broadening terms included according to 
Reardon and Huffaker (Ref. 13) and Reardon et al. (Ref. 14) is 



£ 


/ \ 

/ 273 \ l ' J , / \ 

„ / 273\ 1,1 



. p > { —> + Ki). t 

P i (—j 

(2.46) 


T = 273°K T = 273°K 


The term i is the species being considered, are the species partial 
pressures in atmospheres. The exponent values of m ^ = l/2 and m ^ = 1 
were recommended by General Dynamics/Convair (Ref. 12). The summation, 
j, runs through the number of species in the gas. A representative set of 
the constants needed to calculate the line width with collisional broadening 
for, in this example, the water molecule is presented in Table 2. 
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Table 2 

LINE WIDTH WITH COLLISIONAL BROADENING FOR WATER 


n . . 
i. J 


•5NL '.W 


T = 273 K 


T = 273 K 


n. . = l/2 and n. . = 1 
i,J i.i 


for i y . . becomes 

j ' ij 


j =H 2 ° 

y . . = 
r iJ 

0.09 

= N 2 

=: 

0.09 

(M 

0 

1! 

= 

0.04 

H 2 

=: 

0.05 

= co 2 

z: 

0.12 

= CO 

=: 

0.10 

onating) 


0.44 
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The evaluation of the l/d value to complete the calculation of a, is 
defined by Reardon and Huffaker (Ref. 13) as 



(2.47) 


where 


= 10 


< b „ + % t2 > 


(2.48) 


The and c^ are constant over spectral regions and for water are listed 
as a function of temperature in Ref. 13. 


While 


V = 


0.44 + 0.09 




l/2 


1/2 


C°+ 0.044 (-^) (1 -C°) (2.49) 


where 


c° = -0.1002 + 0.2802 X 10" 3 T - 0. 1089 x 10" 6 T 2 + 0.0291 x 10 _9 T 3 (2.50) 


Values of l/d shall be tabulated as a function of \ and T in Ref. 14. 
These tabulated values of l/d provide an alternate method of obtaining the 
fine structure parameter, a. 

If more than one species is optically active in a given spectral interval, 
Eq. (2.42) is modified and used thusly: 
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called optical depth and the MS indicates multi-species. The type of summation 
indicated in Eq. (2.51) is not obvious, but is what is used. 


The origin of the relationships necessary to determine temperature and 
concentrantions from zone radiometry experiments has now been established. 
These relationships will now be applied to the specific experiments which have 
been performed by Rocketdyne. 
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Section 3 

DATA REDUCTION PROCEDURE 

The solution to the equation of transfer — Eq. (2.29), the Curtis-Godson 
approximation as stated in Eq. (2.37), and specified band model parameters 
may now be used to reduce zone radiometry data. Brewer (Ref. 16) describes 
a computer program to perform such a calculation; unfortunately, he uses a 
distribution function, f , which is not compatible with the reported General 
Dynamics/Convair k's and a's. This may not introduce a significant error, 
but it prevents one from using the reported program directly. 


Rocketdyne chose not to reduce their data in this manner. They approx- 
imated Eq. (2.37) with its grey gas limit Eq. (2.38) and then reduced the radi- 
ometry data, arguing that since a — - oo , X — 0 for CO^ and that experiments 
with water first were optically thin and second that they used water vapor radia- 

A 

tion data which were taken with the same resolution as their spectrometer (Ref. 
15). The first two of these arguments may well be valid, and their validity can 
be determined with analysis. The third is highly improbable because not only 
isothermal property data must be available (which may be) but data for the same 
temperature and compositions as those in the measured plumes must be also. 
However, if the first two arguments are valid, the third is unnecessary. Herget 
(Ref. 15) contends because of these arguments that his studies are not limited 
to optically thin cases. Let us reserve judgment on this contention until some 
subsequent calculations and experimental observations are made. 


Specifically, Eq. (2.38) may be written as 


T k = ( T l) (T 2^ 


n 5 


(3.1) 


th 

where the subscript k represents the k row of zones and h represents the 
number of zones. Eq. (2.29) maybe multiplied by 7r , so that N ^.'s appear on 

22 


LOCKHEED - HUNTSVILLE RESEARCH & ENGINEERING CENTER 



LMSC-HREC D306101 


the RHS where N w ^ is emissive power. In fact, a summation of any of the terms 
listed in Table 1 could have been used. Actually neither intensity nor emissive 
power is measured but the radiation from the solid angle (3 which intersects 
the detector surface. This angle is not measured, but the radiation from an 
internal blackbody source along 3 is. Since any of the terms in Table 1 may be 
calculated from the internal blackbody temperature, the LHS of Eq. (2.29) is 
simply calibrated. Eq. (2.29) is used as stated below. 

+I *b2< T ! T 2> + --- 

= I Abl (e l ) + *Ab2 (T 1 &-T 2 ])+••• 


= 1 


Abl ' 1 


( € i ) + *Ab2 * T 1 + 1 




Abi i T i-1 


Ab3 


* ❖ * ❖ * 

(Tj t 2 . Tj r 2 Tj) 

V * 

sjc # 

T 1 T 2 € 3 


❖ 

t 2 e 3 


(3.2) 


e is defined as 1 - T* , and t ' is defined by Eq. (2.43). Since 1^ is an averaged 
intensity over some small spectral interval and is measurable, it will be called 
radiance. 


Before considering the Rocketdyne experiments in more detail, consider 

❖ 

the following definition of a new T , namely, 


-fn (r**) = 


k h S h 


1 + 


(£ “.A} 2 

n = 1 


n 1/2 


4 y' k S a , 
/ i n n n| 


n = 1 


(3.3) 
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This r h has the property that if it is used in Eqs. (3.1) and (3.2), Eq. (2.42) 
will result. This means that the assumption on optical depth which was used 
by Rocketdyne will be removed. Other means could be used to eliminate this 
assumption, but, as will be confirmed in subsequent discussion, this means 
would require the least amount of revision to the existing Rocketdyne data 
reduction program. An additional benefit is that a convenient check of the 
deviation from optical thinness can also be made with this parameter. 

In summary, the following equation must be solved, either exactly or 
approximately, to reduce zone radiometry data for one component and wave- 
length. Functionality is emphasized for clarity. 


I x KM - £ I xb) X ' T i)' [v, |VA< 0 ,T.. l f-T i )X,Aa,T.l] (3.4) 


Now the specific geometry of the zone radiometry experiments can be 
considered. 

3.1 ONE -DIMENSIONAL ISOTHERMAL FLOWS 

The one -dimensional test situation (Fig. 3) such as the Rocketdyne 
Composite Engine Study (Ref. 17) will be used to demonstrate the procedure 
used to convert measured values of emissive power into temperature and 
composition values. Recalling Fig. 1, the basic geometry of this one- 
dimensional flow contains all of the features shown except that there is a 
single zone. Two radiance readings are made using the zone radiometer. 
One radiance reading is made with the chopper closed giving the radiance 
of only the zone while the other reading made with the chopper open provides 
a radiance value containing the zone radiation and the grey body source radi- 
ation. Using the finite difference form of Eq. (3.2) with these measured 
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radiances, 1^^ (with the chopper closed) and 1 XB (with the chopper open), 
the transmittance and blackbody radiance of the zone can be calculated from: 


[ U= hb (1 - T > 


X XB I \b ^ + *Xb 


(3.5) 

(3.6) 


source 


Note that the temperature of the radiating gas is assumed not be change when 
it is impinged upon by the radiating source. 

The zone radiometry measurement data presented by Rocketdyne are 
in terms of radiance with units of W cm sr jLl . The form of Planck’s 
law applicable for relating the blackbody radiance at a particular wavelength 
to the temperature was used to evaluate the temperature of the zone (see 
Table 1). 



,5 , 

< r (he \ 

,\ 


[| ex P(jxT). 

• V 


(3.7) 


The evaluation of the partial pressure of the radiating species makes 
use of the representation of transmittance. When the species can be treated 
as a continuum radiator the transmittance can be calculated using Eq. (2.45) 

$ 

-In T = k S 

where the S is distance and the k value contains the partial pressure of the 
species as a correction factor to the absorption coefficient at standard con- 
ditions, k Q . Repeating Eq. (2.44) 
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The tabulated k Q value corresponding to the calculated zone temperature 
makes the solution for the partial pressure straightforward. 

Whenever band models must be used to represent the radiation process, 
the transmittance is given by Eq. (2.42). The solution for the partial pressure 
is no longer a simple process. The line width parameter, a, is dependent on 
temperature, the local pressure and partial pressure of all the constituent 
species as described previously. An iterative method is used in which an esti- 
mate of the partial pressure of the radiating spacies is made using the simpler 
continuum radiation transmittance (Eq. (2.45)). With this estimated pressure 
and the calculated temperature, a corresponding "a" value is evaluated using 
Eqs. (2.46-2.50). The first iteration on the partial pressure can then be made 
using Eq. (2.43). The iteration procedure is continued until the desired degree of 

agreement is obtained between succeeding pressure values. 

Summarizing the procedure for evaluating the temperature and partial 
pressure from zone radiometry measurements of a one -dimensional flow: 

1. Obtain the zone radiance and transmittance using Eqs. (3.5) 
and (3.6). 

2. Solve for the zone temperature using Eq. (3.7). 

3. Use Eq. (2.45) for continuum radiation to calculate the 
partial pressure to complete the solution, or 

4. Use an iterative procedure to calculate the partial pressures 
for non-continuum radiation requiring a band model repre- 
sentation of transmittances as follows: 

a. Estimate a partial pressure value using the 
continuum radiation representation for the 
transmittance, Eq. (2.45). 
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b. Evaluate the fine structure parameter, a, 
with Eqs.(2.46) through (2.50) using the 
temperature and the estimated pressure 
value. 

c. Calculate the partial pressure using Eq. (2.43). 

d. Compare the newly calculated partial pressure 
with that used in Step 4b. 

• "Poor" agreement: Repeat from 

Step 4b using new partial pressure. 

• "Good" agreement: Consider solution 

complete. 


3.2 AXISYMMETRIC NONISOTHERMAL FLOWS 

Application of zone radiometry measurement techniques to axisymmetric 
flows uses the same principles as for the one-dimensional situation but the 
reduction of the measured intensities to temperature and partial pressure 
becomes more complex. 

A schematic of the axisymmetric zone layout is given in Fig. 4. 



* 

Zone Number 


Fig. 4 - Definition of Zones and Lines of Sight in Axisymmetric Zone Radiometry 
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The zones consist of concentric circular regions in which the physical 

properties are assumed uniform. The line of sight (L>OS) is defined such 
th th 

that the n LOS passes through the n zone and all zones outside of it. 

The zones are not necessarily of the same size. (In the Rocketdyne zone 
radiometry data reduction program (see Appendix and Ref. 18) the zones are 
nevertheless assumed to all be of the same size.) Whether or not the same 
size zones are used the path lengths within the zones are variable and depend- 
ent upon the location of the zone within the axisymmetric array. The path 
length is calculated using geometric relationships. For example the path 
length, £ , in the fourth zone on the third line of sight is (see Fig. 5) calcu- 
lated as 

i = y Rj - r j 


V 


R 2 " R2 
C 1 


(3.8) 



Zone Number 


Fig. 5 - Path Length Definition in an Axisymmetric Zone Radiometry Array 
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From Fig. 4 it can be seen that the signal received by the detectors 
along a respective LOS has in general passed through an inhomogeneous 
region. The radiances measured for the LOS can be used to calculate the 
temperatures and partial pressure within the zones using the following 
procedure. 

Again the two radiance readings (one with the chopper and one without 
it) are made along each line of sight. The zones are maintained at sufficiently 
small sizes that the line of sight through the concentric zones can be approxi- 
mated by one-dimensional slabs as in Fig. 6. 



Zone Number 


Fig. 6 - One -Dimensional Approximation of the Lines of Sight for Axisymmetric 
Zone Radiometry 


The radiance values measured for these lines of sight can be mathe- 
matically represented as they were in Eqs. (3.5) and (3.6) for the single one- 
dimensional case. 
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r XAj ^bi, j € i ’ J T i-l,j 


(3.9) 


i = l 


^XPj ]^*Abi, j e. , T i-l,j * *Ab Source T j 
i=l 1,J 


(3.10) 


where j is the line of sight under study and the i is summed over all 
the zones to n. Subtracting Eq. (3.9) from Eq. (3.10) provides n relationships 
for 7. which is the mean transmittance of the entire 4 th line of sight. The other 
n equations needed to solve for the 2h unknowns, and P. come from Eq. (3.9) 
The representation of the transmittance now becomes the prime question. 
Rocketdyne uses the grey gas limit for the transmittance, Eq. (3.1). The n 
values of Tj can then be expanded as 


n * * * $ 

T ] = n r = (r .) (r 2 .) ... (T ) 

J , i = 1 1 > J 1« J J n » J 


(3.11) 


The 2n equations consisting of Eqs.(3.9) and (3.11) are solved using matrix 
algebra. Since Rocketdyne has automated the solution procedure in a data 
reduction program, an iterative process is used to evaluate the unknowns. 

Summarizing the procedure for evaluating the temperatures and partial 
pressures in the n zones of an axisymmetric flow using the zone radiometry 
measurements is: 


1. Construct equations for the measured line of sight radiance 
using Eq. (3.9). 

2. Construct equations for the mean transmittance through an 
entire line of sight using Eq. (3.11). 

3. Place the transmittance represented by Eq. (3.1) in matrix 
form. (Using the measured mean line -of-sight transmittance 
values, T- , the matrix can be solved for the zone transmit- 
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4. Replace the transmittance s in the equations constructed 
in Step 1 with the calculated zone transmittances. (The 
matrix representing the mean measured radiance is then 
ready for solution for the zonal blackbody radiance func- 
tions.) 

5. Solve for the temperatures in the zones using Eq. (3.7). 

(For continuum radiators this completes the solution 
procedure since the zone partial pressures can be eval- 
uated using the calculated k. 's, tables of k versus 
temperature and Eq. (2.44).) 1 u 

6. Use an iteration procedure (for noncontinuum radiation, 
requiring a band model representation of the transmittance) 
to solve for the temperatures and partial pressures in the 
zones as follows: 

a. Use Eqs. (2.46) and (2.50) to evaluate the fine structure 
parameter, a, using the zone temperatures and the 
zone partial pressures from Step 5. 

b. Reevaluate the zone transmittance values using Eq. (2.43) 
The modification to the Rocketdyne program to define 
the transmittance of the zone with Eq. (3.3) would elim- 
inate the grey gas assumption inherent in Eq. (2.43) and 
make the solution valid for all optical thicknesses. 

c. Return to Step 4 and repeat Steps 4 and 5. 

d. Compare the newly calculated partial pressures and 
temperatures of the zones with those obtained previously 
in Step 5. 

• "Poor" agreement: Repeat from Step 6 using new 
partial pressures and temperatures for the zones. 

• "Good" agreement: Consider solution complete. 


3.3 ROCKETDYNE ZONE RADIOMETER DATA REDUCTION PROGRAM 


The Rocketdyne automated data reduction program is listed in the Appen 
dix. An input guide and flow chart of the program are also presented. To aid 
potential users of the data reduction program, a sample case is given along 
with sample input and output. 
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Section 4 

EXAMPLE PROBLEMS 

To illustrate the calculation techniques discussed in this report, several 
example problems will be solved. The first is one typical of an axisymmetric 
alcohol -burning Atlas vernier engine; the second represents a planar, hydrogen 
oxygen diffusion flame, i.e., the composite engine experiment. 

• Problem 1 — LOX- Alcohol Engine 

Due to the behavior of gaseous radiation properties, it is desirable to 
choose example problems in which the temperature and composition are spec- 
ified. Consider first the following thermal path: 


T (°K) 


Detector 


1500 

2000 

2500 

2000 

1500 

— 

3500° source 

Zone: 

1 

2 

3 

4 

5 




Zones 1, 2, 4 and 5 are 2 cm long; zone 3 is 4 cm long. 

Partial pressure, CO^- 0.27 atm 

H^O: 0.58 atm 

Wave lengths of measurement: 4.45 (/i) or 2247 (l/cm) 

2.49 (/i) or 4016 (l/cm) 


Solution Procedure 
1 . Evaluate 


k = k 


273°K 
T°K 


r P. atm' 

j. 

1 atm 
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k Q (l/cm) 

T (°K) 

oo (l/ cm) 

Ref. 

Species: 

1500 

1 

2000 

2500 

2247 

4016 

(12), p. 70 
(12), p. 73 

h 2 o 

1.22 x 10’ 2 
1.47 x 10 _1 

2.33 x 10" 2 
1.43 x 10 _1 

3.05 x 10 -2 
1.50 x 10" 1 

C °2 

10.99 

(Not measv 

i 

13.30 

ired but assui 
! 1 

13.55 
med zero) 

2247 

4016 

(10), p. 104 
(12), p. 33 


T (°K) 

T/273 (°K) 

[t/2.73 (°K) 

1/2 


T 3 

mmm 

5.50 

2.35 


2.25 x 10 6 

3.38 x 10 9 

mSm 

7.34 

2.73 


4.0 x 10 6 

8.0 x 10 9 

2500 

9.15 

3.03 


6.25 x 10 6 

15.63 x 10 9 


k (l/ cm) 

T (°K) 

00 (l/ cm) 

Species: 

1500 

2000 

2500 


H,0 

1.29 x 10 -3 

1.84 x 10' 3 

1.94 x 10' 3 

2247 

L 

1.55 x 10 -2 

1.12 x 10 -2 

0.95 x 10 -2 

4016 

C °2 

0.540 

0.485 

0.400 

2247 


2. Evaluate fine structure parameters. 


a. line density, (l/d) 

(i) C0 2 , (l/d) = (l/DLR) (cm) at 00 = 2247 (l/cm) 


(l/DLR) 

T (°K) 

181.1 

1500 

356 

2000 

510 

2500 


Ref. 11, p. 59 
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(ii) H z O, (l/d) = a */y* Ref. 13, p. 194 

£ 

Evaluation of a is a rather nebulous operation, but Rocketdyne's pro- 
grammed values will be presented. Reference 14 is supposed to have tabu- 
lated values of (l/d) when it is published. 



To complete the (l/d) calculation for H 2 O, define y = | 0.44 (T q /T) + 
0.09 (T q /T) 1//2 ] C* + 0.044 (T q /T) 1 ^ (1 - C*) Ref. 13, p. 194 

where 

C* = -0.1002 + 0.2802 x 10“ 3 T - 0.1089 x 10" 6 T 2 
+ 0.0291 x 10" 9 T 3 . 
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Evaluate C and y for the temperatures in the zones. 

For T = 1 500°K 

C* = -0.1002 +0.423 

-0.245 +0.098 

-0.345 +0.521 = +0.176 

T = 2000°K 

C* = -0.100 
-0.436 
-0.536 

T = 2500°K 

C* = -0.100 +0.701 

-0.682 +0.456 

-0.782 +1.157 = +0.377 

y* = [o.44 (T /T) + 0.09 (T o /T) 1/,Z ] C* + 0.044 (T q / T) 1 ^ 2 (1 - C*) 
Use the C values to calculate the y values. 

y*50o = [0.080 + 0.0383] (0. 176) + (0.0187) (0.824) 

= (0.1183) (0.176) + 0.0154 

= 0.0209 + 0.0154 = 0.0363 

y*ooo = [0*0593 + 0.033] (0.257) + (0.0161) (0.743) 

= 0.0237 + 0.0120 = 0.0357 

72500 = [0.0481 + 0.0297 [0.377] + (0.0145) (0.623) 

= 0.0292 + 0.0091 = 0.0393 
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The (l/d) values for H^O at the temperatures and in the wave lengths of 
interest are: 


(l/d) (cm) 

T (°K) 

i (l/cm) 

2.78 

1500 

2247 

5.49 

2000 


11.78 

2500 


1.68 

1500 

4016 

3.52 

2000 


8.15 

2500 

1 

1 


b. Obtain the collision half widths, y^'s, to complete the calculation of 
the fine structure parameters, a's. 



j = all species. 


For CO z at 1500°K, Ref. 14 


y = (0.07) (0.58)/2.35 + (0.09) (0.27)/2.35 + (0.01) (0.27)/5.50 

C °2 

= 0.0173 + 0.0103 + 0.0005 = 0.0281 (l/cm) 


At 2000°K, 

y = (0.07) (0.58)/2.73 + (0.09) (0.27)/2.73 + (0.01) (0.27)/7.34 

CO, 

= 0.0149 + 0.0089 + 0.0004 = 0.0242 
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At 2500°K 


v 



(0.07)(0.58) ,. + _j0 J) 9I (0 , : 27) , + (0 . 0 1 > (0.27)/9. 1 5 


(0.0406 + 0.0243) 
3.03 


0.0003 = 0.0217 


For H z O: Ref. 14 


y = l °-°9) (0-58) + (0 12) (0.27J + ( 0 .44) (0. 58)/(T/273) 

c h 2 o 

0.0845 0.255 

(T/273) 1 ^ 2 (T/273) 

At 1 500°K 

y = 0.0360 + 0.0465 = 0.0825 

c h 2 o 

At 2000°K 

y = 0.0310 + 0.0347 = 0.0657 

c H 2 o 

At 2500°K 

y = 0.0279 + 0.0279 = 0.0558 

c h 2 ° 


The radiation parameters for this sample problem are summarized on 
the following page. 
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a = (r c )(l/d) 

T(°K) 

Species 

0) (l/ cm) 

k (l/ cm) 

5.09 

mm 

CO- 

2247 

0.540 

8.61 

■SEE 1 

L 


0.485 

11.10 

: smm 



0.400 

0.229 

0.358 

0.658 

1500 

2000 

2500 

H 2 o 

2247 

1.29 x 10"? 
1 .84 x 10 ~ 
1.94 x 10~ 

0.139 

0.230 

0.455 

1500 1 

2000 
2500 

H 2° 

4016 

1.55 x 10_? 
1 . 1 2 x 10 
0.95 x 10"* 


The intensities measured by the detector are the blackbody radiation of 
the species within the zones attenuated by the zones between that particular 
zone and the detector. 


The blackbody intensity for the zonal temperatures and the wave numbers 
of interest are < 


hb (Ret - 19 > 

T 

i 

0.894Z (W/cm^- n-sr) 

1500 (°K) 

2247 (l/cm) 

1.691 

2000 


2.580 

2500 


4.493 

3 500 


2.687 

1500 

4016 

7.271 

2000 


13.56 

2500 


29.19 

3500 



The attenuation (or in the converse sense, the transmittance, T) of the 
radiation of the zones to the detector is calculated for each of the optical paths. 
Starting from the detector, the first transmittance, T ,, involves only zone one. 
The calculation of the second transmittance, T 2 » includes the first and second 
zones. These calculations continue until all the transmittances are determined. 
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For CC >2 radiation: 

Kj Lj = (0.540 cm ^)(2cm) = 1.080 



Tj = 0.350 

K 2 L 2 = (0.485) (2) = 0.970 



T 2 = 0.144 

K 3 L 3 = (0.400) (4) = 1.60 
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- £nT 3 


2.05 + 1.60 


T 72 


1 + 


(3-65)' 


[4(13.85 + (1 1 . 1 ) (1 . 6)]y 


3.65 

(1.052) 


3.48 


T 3 = 0.0310 


K 4 L 4 = K 2 L 2 =0 - 970 


-£n r 4 



3.65 + 0.97 
(4.62) 2 

1.65 + ( 8. 6 1 ) (0 . 



TJl 


4. 62 
1.06 


4.35 


r 4 = 0.0130 


K C L C = K.L. = 1.080 

D O 11 


-£nr 5 = 



4.62 + 1.08 


T7T 


(5-70)^ 


4 [40.00 + (5.09K1.08)] 


5.70 

1.09 


5.25 


T c = 0.0053 

D 


Use the zonal blackbody radiation and transmittance values to evaluate 
the measured intensity. 




Abl 


U-rp 


+ \bZ * T 1 


T 2> 


+ I Ab3 (r 2' T 3 ) + I Ab4 (r 3- 




+ ^ 5^5 ' 
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I x = / 0.894(1 - 0.350) = (0.894) (0.650) = 0.580 ^ 

1.691 (0.350 - 0.144) = (1.691) (0.206) = 0.348 

2.580(0.144- 0.0310) = (2. 580) (0. 1 13) = 0.292 

1.691 (0.0310 - 0.0130) = (1.691) (0.018) = 0.031 
0.894(0.0130 - 0.0053) = (0.894) (0.008) = 0.007 ) 


Ix = 1.258 w/cm -/tl-sr 


If a 3500°K blackbody source has also been transmitting, 

I xws = 1.258 + (0.0053)(4.493) 

= 1.272 w/cm^-jU-sr 

Neglecting non-grey effects introduces errors of 5 to 10% in intensity; these 
become a smaller precentage when they are converted to temperature. Such 
errors do not become unreasonably increased when using the data reduction 
programs — as evidenced by Appendix D of Ref. 18. 

For H^O radiation at 2247 (l/cm), the transmittances are 


KjLj = 2.58 x 10 


■fnTj = 


2.58 x 10 


[ . 2.58 x 10~ 3 \ 

^ + 4(0.229) f 


172 


2.58 x 10" 3 „ ,„-3 

= 2.5o x 10 

(1. 0028/' C 


rsJ . 

7 1 = 1 


K 2 L 2 


= 3.68 x 10 


-fn T 0 = 


6.26 x 10 


-3 


1 + 


(6.26 x 10~ 3 ) 


772 


6.26x 10 3 _, lrt -3 
1.0026 6 - 26x 10 


4 [o. 59 x 10" 3 + 1.32 x 10’ 3 L 
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T 2 = 1 


K 3 L 3 = 7.76 x 10 


-3 


-£nr 3 = 


14.02 x 10 


-3 


1 + 


-2 

( 1 .40 x 10 ) 




4[l.9lx 10 " 3 + 5. 1 1 x 10 -:5 ] 


-3- 


l,40x 10 
1.005 


-2 


t 3 = 0.9861 


K 4 L 4 = K 2 L 2 = 3.68 x 10 


-3 


-*nr 4 = 


17.70 x 10 


-3 


1 + 


(1.770 x 10 ~ 2 ) 2 


172 


1.77x 10 
1.005 


-2 


4 [7.02xl0 _3 + 1.32xio -3 ] 


T 4 = 0.9825 


K 5 L 5 = K J L 1 = 2.58x 10 




20.28 x 10 


-3 


1 + 


400 x 10 


-6 


172 


2.03 x 10 
1.005 


-2 


4 [8.34x 10 ~ 3 + 0. 59 x 1 0 - " 5 ] 


-3- 


T 5 = 0.9799 


X Xb2 (T 1 ' t 2 ) = 0 

^3 (r 2 - t 3 ) = (1 - 0.9861) = 0.0139 
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r Xb4 < T 3 " T 4> = (0 - 9861 ' 0.9825) = 0.0036 
hb5 (T 4 ' T 5 ) = (0 - 9825 ' 0.9799) = 0.0026 


L = 0.892 (0) + 

A 1.691(0) + 

2.580 (0.0139) + 
1.691 (0.0036) + 
0.8921 (0.0026) 


0 

0 

0.0358 

0.0061 

0.0023 


l x = 0.0442 


w/cm - li- sr 


The intensity when the 3500°K blackbody source is also transmitting is 


X XWS = °* 0442 + (0.9799)(4.493) 

= 4.434 w/cm^-jLi-sr 

In this example, water vapor radiation is very optically thin. 

Total Radiation at 2247 (l/cm) 

Since both CO^ and H^O are optically active at 2247 ( 1 / cm), both contribute 
to the total radiation flux. To account for multi-species emission, the natural 
log of the transmittance of each is calculated and then all such logs are summed. 


-£n (t, ) 

MS 

-fn(T ? ) 

MS 

-jMTo) 

MS 

-*n(r 4 ) 

MS 

-ln(r,) 

MS 


1.050 + 0.0026 = 1.053 
1.940 + 0.0062 = 1.946 
3.48 + 0.014 = 3.49 
4.35 + 0.018 = 4.37 
5.25 + 0.020 = 5.27 
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(T h ) 

MS 

<T - ‘ ^MS 


0.348 

0.652 

In this case the 

0.143 

0.205 

results are identical 

0.0305 

0.112 

to those for CC> 2 

0.0126 

0.0179 

alone. 

0.0051 

0.0075 



Since the radiance from the water vapor is only 4% of that from C0 2 » 
Rocketdyne assumed it negligible (consistent with the exact calculation). 

Repeat the transmittance and intensity calculations for H^O radiation 
at 4016 (l/cm) 


= (1.55 x 10' 2 )(2) = 3.10 x 10 -2 


a A = 0.139 


-fnTj = 


3.10 x 10 


1 + 


3.10 x 10 
4(0.139) 


772 


3.10 x 10 = 3>02 x 1q -2 


1.025 


Tj = exp (-0.031) = 0.9704 
K 2 L 2 = 2.24 x 10" 2 


-£nr 2 = 


5.34 x 10 


1 + 


28.4 x 10 


-4 


772 


5.34xl0' 2 , ..-2 

OJ3— = 5.13x10 


4[4.31x 10 -3 + 5. 15x 10” 3 }/ 


r 2 = 0.950 


K 3 L 3 = 3.80 x 10 


-2 


-£nr 3 = 


9. 14x 10 


-2 


1 + 


83 x 10 


771 


= 8 , 80xl0 - 2 


4 [9.46x 10" 3 + 1.73x 10“^] 


- 2 - 


45 

LOCKHEED - HUNTSVILLE RESEARCH & ENGINEERING CENTER 




LMSC-HREC D306101 


t 3 = 0.916 


K.L. = 2.24 x 10 
4 4 


-2 


-in r 4 = 


11.38 x 10 


-2 


1 + 


130 x 10 


-4 


172 


1.138x 10 _ 1>08x 1Q -1 


1.05 


4[26.76x 10 -3 + 5.15 x 10” 3 ]j 


r 4 = 0.898 


K 5 L 5 = 3.10 x 10 


-£nT 5 = 


14.48 x 10 


-2 


1 + 


210 x 10 


-4 


4 [3 1.9 lx 10" 3 + 4.3 lx 10"^ J 


_ 1 .448 x 10' 1 , „ r , „-l 

172 TT07 = 1.35x10 


“3n 


T 5 = 0.874 


L Xb 


X Xbl 

(1- 

T l> = 

(2.687) 

(0.0296) = 

0.080 

X Xb2 

(T 1' 

" T 2> = 

(7.271) 

(0.0204) = 

0.150 

X Xb3 

(t 2' 

■r 3 ) = 

(13.56) 

(0.034) = 

0.460 

X Xb4 

(t 3 - 

-t 4 ) = 

(7.271) 

(0.018) = 

0.131 

X Xb5 

(t 4‘ 

" T 5> = 

(2.687) 

(0.024) = 

0.065 


I x = 0.886 w/cm 2 - jLt-sr 

I> = 0.886 + (0.874)(29.19) = 0.886 + 25.50 = 26.39 W/cm 2 -jU-sr 
A WS 
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• Problem 2 — Composite Engine 

This is an example of single-zone radiation. Let the path length be the 
distance between the side walls less the initial jet widths of coolant gases, 
4.46 — 0.80 = 3.66 in. or 9.30 cm. Static pressure 15.3 psia. Mass fractions: 
0.9438, water and 0.0545, hydrogen. 

The mole fraction of water is 0.656, giving a partial pressure of 0.683 
atm. The wave number of interest is 4016 (l/cm). 


k[l/cm) 

T (°K) 

kL 

(1 + kL/ 4a) 1 / 2 

r 

1.825 x 10" 2 

1500 

0.1695 

1.145 

0.8630 

1.33 x io"; 

2000 

0.1235 

1.065 

0.8910 

1.12 x 10" ^ 

2500 

0.1040 

1.030 

0.9040 


The measured intensity for the zone can be represented by 1^ = 
In particular for each zone 


(0.137) (2.687) = 

0.368 

at 

1 500°K 

(0.109) (7.271) = 

0.792 

at 

2000°K 

(0.096) (13.56) = 

1.310 

at 

2500°K 


When a 3 500°K blackbody source is also radiating, the measured intensity 
for each zone is calculated as 

x xws = 0,368 + (0,863)(29,19) = 0,368 + 25,05 = 25,418 at 1500 ° K 

= 0.792 + (0.89l)(29. 19) = 0.792 + 25.95 = 26.74 at 2000°K 

= 1.310 + (0.904){29.19) = 1.310 + 26.20 = 27.51 at 2500°K . 
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The experiment corresponding to this calculation would indicate: 


hMS - 27 51 
I X = 1.310 


therefore, 


27.51 - 1.310 
29.19 


26.20 

29.19 


0.900 


1 - T * 0.100 

X Xb = ~ 5T6 ' 6 = 13 - 10 » T = 2470°K = 2500°K • 


Not using band models would introduce intensity errors of up to 15% in 
the variable range presented here; temperature errors would be somewhat 
less. If € is much less than 0.1, serious errors would be introduced into the 
temperature determination. 


The assumption of a grey gas may be used for CO^ and of an optically 
thin gas may be used for P^O, in the examples presented, without introducing 
excessive errors. Such errors could be removed by using a more complete 
data analysis program. The theoretical radiation analysis presented in Section 
2 should provide a very adequate basis for an accurate data analysis calculation 
in ranges of experiments for which the illustrative examples are typical, i.e., 
no improvement to the Curtis-Godson approximation is necessary. 
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Section 5 

CONCLUSIONS AND RECOMMENDATIONS 

This report has demonstrated that sufficient radiation property data 
exist to study zone radiometry of CO^ and H^O. Such data also exist for 
soot and CO. 

Data reduction procedures currently used are adequate for the studies 
which Rocketdyne has performed. These work because CO^ is a grey gas 
and H^O is optically thin in their experiments. Sample problems show this 
behavior. More accurate data reduction schemes could be devised, but this 
would not substantially improve existing data. However, to remove that 
criticism such procedures should be developed. 

The only apparent reason for experiments of the alcohol-LOX vernier 
engine type not yielding accurate temperature and partial pressure data is 
lack of axial symmetry. The two-dimensional mixing study is probably so 
optically thin, that accurate transmittances would be very difficult to deter- 
mine. In any event, all future studies should be preceded by an error analys 
of expected data. 
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Appendix 

The data reduction procedure for axisymmetric zone radiometry read- 

* 

ings was automated by Rocketdyne in a computer program called the 
Axisymmetric Zone Radiometer Data Reduction Program. A description 
of the program is given here by outlining the operations occurring in the sub- 
routines, a detailed flow chart of the procedure and a listing of the program. 
An input guide is given, and also a sample input and output are included to aid 
potential users of the program. 

PROGRAM SUBROUTINES 

MAIN Subroutine 

The Axisymmetric Zone Radiometer Data Reduction Program main 
driver is responsible for reading the program input and calculating the local 
temperature and species concentrations of water or carbon dioxide. 


The main program does the following (in sequence): 

1. Evaluates the path lengths 

2. Reads data 

3. Solves for the product of zonal absorption coefficient and zonal 
partial pressure of radiating species, (kp)^ 

4. Constructs attenuation matrix to solve for zonal blackbody 
radiance 

5. Solves for the zonal blackbody radiance 

6. Uses Planck's distribution law to solve for zonal temperatures 

7. If the radiating species are considered to be continuous radiators, 
the solution is complete. The local temperatures along with the 
product of the absorption coefficient and partial pressures for 
each zone is output. 


North American Rockwell Corp., "A Compendium of Zone Radiometry 
Measurements of Exhaust Plumes," R-8140, Rocketdyne Div., Canoga 
Park, Calif., 25 February 1970. 
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8. If the radiating species require the use of spectral averaged 
data (band models), the zonal temperatures and partial pres- 
sures are iterated. 

9. Evaluate the fine structure parameter, a, for each zone using 
the zonal temperatures and partial pressures. 

10. Calculate new absorption coefficients for each zone. 

11. Repeat calculation steps from Step 4 until successive values of 
the temperatures and partial pressures are within a preset limit. 

QUAD Subroutine 

This subroutine contains three entries: QUADO, QUAD1 and QUAD2. 

The: QUADO entry generates the appropriate weighting function for quadratic 
distribution of properties. The QUAD1 entry sums the product of path lengths 
and weighting factors to obtain the coefficient for the average (kp) values. The 
QUAD2 entry sets up the average kp values using the path lengths and the 
weighting factors. 

SLIT Subroutine 

A routine to correct for spatial resolution in deflection data. The 
correction technique is from the University of Tennessee described in AF 
CRL 465, pages 59-62. These corrections are usually small and occur 
at the edge of the plume. Rocketdyne has modified the experimental pro- 
cedure such that radiance and transmission (smoothed) data are available 
rather than deflection data making it unnecessary to apply the slit corrections 
in this subroutine. For completeness the capability to read in deflection data 
and correct it for spatial resolution has been left in the program. 

ISIMEQ Function Subroutine 

This subprogram solves a set of simultaneous linear equations with up 
to 30 variables. This is a standard matrix solution subprogram. Throughout 
the zone radiometry data reduction program, this subroutine is used to solve 
for the variable of interest in each zone and then returns the answers as a 
column matrix in column one. 
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PROGRAM FLOW CHART 
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Path Length Evaluation 


Input Data 
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Call QUAD 2 to get 

matrix of average 

values of k. P, f. . 

1 k. 1 , l 

i J 


fr 
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Solution Complete Write Out 


1. Title 

2. Zone 

3. Average k Product 
at Each Zone 

4. Emissivity for Each 
Zone 

5. Radiation Energy at 
Each Station 

6. Temperature at Each 
Station 

7. Zone Width in cm 
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INPUT GUIDE TO THE DATA REDUCTION PROGRAM 
FOR AXISYMMETRIC ZONE RADIOMETRY 


Card 

Col. 

Format 

Description 

I 

1-72 

18A4 

Title card 

2 

1-12 

112 

Number of zones 


13-24 

F12.8 

Zone width, cm 


25-36 

F12.8 

Correction to plume radiance due to window 




absorption 


37-48 

F12.8 

Blackbody radiance, W/cm^-ji-sr 


49-51 

13 

Property variation 




= 0 Quadratic property variation 




= 1 Step function; constant zonal properties 


52-54 

13 

Slit function correction control parameter 




= 1 Apply a new function 




= 0 Apply previously used function 




= -l Do not correct for slit function 


55-60 

16 

Data type 




= 0 or 1 For spectral data 


| 


= 2 For spectrally averaged data 




(Positive values imply smoothed input — 




zero or negative values imply deflection 




values input) 

2 

61-72 

F 12.8 

Wavelength in microns 

3a 



Deflection Data (or) Smoothed Data 




For Deflection Data 


1-12 

F12.8 ' 




13-24 

F12.8 


Blackbody Data; 6 values to a card until 


— 

F12.8 


all zone data are input 


61-72 

F12.8 



3b 

1-12 

F12.8 : 




13-24 

F12.8 


Flame Data; 6 values to a card until 


j 

F12.8 


all zone data are input 


61-72 

F12.8 

> 
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Card 

Col. 

• 

Format 

Description 

3c 


F12.8 




1 

F 12.8 


Blackbody Prime Data; 6 values to a card 


B 

F 12,8 


until all zone data are input 


61-72 

F12.8 



3d 

1-12 

F12.8 




13-24 

F12.8 


Flame Prime Data; 6 values to a card 


— 

F12.8 

> 

until all zone data are input 


61-72 

F12.8 



3e 

1-12 

F 12.8 ' 




13-24 

F12.8 


Greybody Alone Data; 6 values to a card 


— 

F12.8 


until all zone data are input 


61-72 

F 12.8 



3f 

1-12 

F 1 2 . 8 




13-24 

F 12.8 


Greybody through Flame Data; 6 values 




F12.8 

► 

to a card until all zone data are input 


61-72 

F 12.8 



' 3g 

1-12 

F 12.8 




13-24 

F12.8 


Greybody Alone Prime Data; 6 values to 


— 

F12.8 

► 

a card until all zone data are input 


61-72 

F12.8 

1 

1 


3h 

1-12 

F 1 2 . 8 ' 

! 



13-24 

F 12.8 

i 

Greybody through Flame Prime Data; 


— 

F 12.8 


6 values to a card until all zone data 


61-72 

F12.8 


are input 

3 a 

1-12 

F12.8 * 


For Smoothed Data 


13-24 

F12.8 


Line of Sight Spectral Radiance Data in 


— 

F12.8 

► 

2 | 
W/cm -micron-ster ; 6 values to a card 

until all zone data are input 


61-72 

F12.8 


3 b 

1-12 

F12.8 




13-24 

F 12.8 


Line of Sight Transmittance Data which is 


— 

F12.8 


dimensionless; 6 values to a card until all 
zone data are input 


61-72 , 

F 12.8 

j 
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SAMPLE INPUT TO THE DATA REDUCTION PROGRAM FOR AXISYMMETRIC 
ZONE RADIOMETRY 


SAMPLE CASE FOR ZONE RADIOMETRY C02 RADIATING 


16 0.48 


• 229 

1.219 

1.196 

.953 

.823 

.647 

• 09 

.033 

.004 

• 18 

.182 

• 1 84 

• 263 

.31 

.372 

• 784 

.92 

.99 



1 -1 

1 4 • 45 

1.163 

1.114 

1 .05 

• 442 

CD 

(VJ 

• 

• 167 

• 00 1 



. 19 

.203 

• 227 

• 453 

.849 

• 66 1 

• 999 
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K 

T<K) 

sample 

NF (F ) 

1 

. 1 BOO 

1 .2290 

2 

. 1870 

1 .2190 

3 

.1890 

I .1960 

9 

.... . 1900 

1 .1630 

5 

• 2030 

1.1190 

6 

.2270 

1 .0500 

7 

. 2630 

.9530 

8 

.3100 

.8230 

9 

.3770 

.6970 

I 0 

.9530 

. .9920 

11 

.5990 

.2800 

12 

.6610 

. 1 670 

13 

.7890 

.0900 

1*» 

♦ 9200 

.0330 

15 

.9900 

.0090 

1& 

.9990 

.0010 


CASE FOR 7 0 N F RAD10MFTRY 


C 02 radiating 


SAMPLE CASE FOR 20ME RADI OMFTRY C02 RADIATING 


radial properties 


F 

<P(F ) 

EPS ( K ) 

RADIK) 

T I F > 

1 

. 1 800 

.0778 

1 .966fl 

2159. 

2 

. 1 793 

.0759 

1.9930 

2 196. 

3 

.18 2 3 

.0787 

1.883o 

2088. 

9 

.1895 

• O818 

1.7679 

2095. 

5 

.1922 

.0879 

1 .699*f 

20C 2. 

6 

.185| 

.0799 

1.6733 

I 990. 

7 

. 1 707 

.07*10 

1 .6299 

196! . 

8 

.1551 

.0679 

1 .5330 

1906. 

9 

. 1 366 

.0596 

1 .3270 

1781. 

10 

.1190 

.0500 

.9889 

1569. 

1 1 

.0917 

.0909 

.7069 

1 Y67. 

1 2 

.O687 

.0305 

.5231 

1 229. 

13 

.0958 

• 02p9 

.9173 

1 ! 33. 

19 

.0169 

.0076 

.9121 

1128. 

15 

.0020 

.0009 

.3736 

1093. 

16 

.0007 

.0001 

1 .0000 

1572. 


70NE 

WIDTH s .950 CM 

USED TO 

COMPUTE EPS 


STEP 

function property 

distribution assumed 
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